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INTRODUCTION 


In  both  the  compressor  and  turbine  the  interaction  of  the  rotor  with  the 
stator  is  critical  to  the  function  of  the  stage.  In  the  rotor  work  is  either 
performed  on  or  extracted  from  the  gas  and  in  the  stator  a  rotational 
component  either  removed  or  imparted  to  the  flow.  In  view  of  the  relative 
motion  of  the  rotor  to  the  stator,  the  basic  interaction  is  unsteady  in 
nature.  In  recent  years  there  has  been  conflicting  evidence  on  the  effect  of 
spacing  on  stage  performance.  Here  spacing  refers  to  the  distance  between 
the  trailing  edge  and  the  leading  edge,  or  vice  versa.  Obviously  there  are 
profound  beneficial  structural  effects  which  follow  from  having  a  short 
spacing  distance.  The  overall  machine  is  shorter,  bearings  have  lower  steady 
loads  and  main  shafts  can  be  shorter  and  stronger  for  a  given  weight.  As 
might  be  expected,  moving  the  rotor  and  stator  closer  together  can  result  in 
a  very  strong,  complex  aerodynamic  interaction.  Moving  the  rotor  and  stator 
further  apart  might  be  expected  to  reduce  the  transient  aerodynamic  loads 
involved,  a  beneficial  effect  obtained  at  a  very  significant  cost.  However, 
the  complexity  of  the  rotor-stator  interaction  has  resulted  in  conflicting 
experimental  evidence  on  the  effect  of  spacing.  The  complexity  of  the 
interaction  is  significantly  greater  in  modern  highly  loaded  or  high  work 
stages  operating  in  the  transonic  regime.  The  indications  are  that  the 
complex  interaction  between  the  stator  and  the  rotor  is  not  as  well 
understood  as  would  be  desired.  An  improved  understanding  could  clearly  have 
a  pronounced  beneficial  effect  on  stage  performance  and  durability. 

In  a  series  of  AFOSR-sponsored  programs,  Detroit  Diesel  Allison  has 
studied  experimentally  various  aspects  of  the  unsteady  flow  in  a  compressor 
(Ref.  I, 2, 3, A)  and  made  comparisons  with  inviscid  theories  of  varying  degrees 
of  complexity.  While  very  encouraging,  much  remains  to  be  done,  particularly 
since  inviscid  theories  cannot  treat  vorticity  production  effects.  Another 
experimental  effort  focusing  upon  the  interaction  problem  is  that  of  Dring, 
Joslyn,  Hardin  and  Wagner  (Ref.  5)  who  investigated  the  effect  of  gap  on  and 
in  the  interaction  process.  These  experiments  clearly  show  the  importance  of 
understanding  the  effect  of  rotor-stator  gap  spacing  and  an  understanding  of 
the  basic  fluid  mechanics  involved  in  this  problem  would  be  a  major  aid  to 
the  research  and  development  engineer. 


In  developing  an  analysis  to  predict  the  stator-rotor  flow  field,  it  is 
clear  that  viscous  effects  must  be  included.  The  work  under  the  present 
effort  attacks  this  problem  by  applying  a  time-dependent  Navier-Stokes 
analysis  to  the  stator-rotor  flow  field.  The  analysis  has  been  successfully 
used  for  calculating  the  flow  about  a  variety  of  cascade  configurations 
(e.g.  Refs.  6-10).  Under  the  present  effort,  the  transient  accuracy  of  the 
Navier-Stokes  code  would  be  assessed,  possible  radiative  boundary  condition 
formulations  included,  a  stator-rotor  grid  developed  and  a  demonstration  case 
run. 

ACTIVITY  DURING  THE  PAST  YEAR 
Transient  Accuracy 

Although  the  Navier-Stokes  code  had  been  well  exercised,  most  of  the 
effort  had  focused  upon  calculation  for  flow  in  steady  configurations.  In 
order  to  apply  the  code  to  the  stator-rotor  problem,  a  time-dependent 
application  is  required.  Therefore,  the  first  task  in  the  program  was  to 
obtain  an  assessment  of  the  time-dependent  accuracy  of  the  code.  In  this 
regard,  two  calculations  were  considered.  The  first  was  done  under  contract 
to  NASA-Langley  Research  Center  as  part  of  a  dynamic  stall  analysis  and  since 
it  is  relevant  to  the  present  effort,  the  results  will  be  presented  here. 

The  full  details  of  the  case  can  be  found  in  Ref.  11.  The  case  considered  is 
that  of  an  NACA  0012  airfoil  oscillating  sinusoidally  in  pitch  in  a  stream 
with  a  Reynolds  number  based  on  chord  of  2.08  x  106  and  a  Mach  number  of 
0.30.  The  airfoil  oscillated  between  4°  and  20°  with  a  dimensionless 
frequency  of  0.125.  The  case  corresponds  to  Data  RUN  51.005  of  the  data  of 
St.  Hilaire  and  Carta  (Refs.  12  and  13).  Although  this  calculation 
represents  an  Isolated  airofll  rather  than  a  cascade,  it  was  made  with  the 
cascade  deck  with  boundary  conditions  appropriate  to  the  isolated  airfoil 
applied. 

Comparisons  between  calculated  and  measured  surface  pressure 
coefficients  are  shown  in  Figs.  1-8.  three  comparisons  during  the  upstroke 
are  shown  in  Figs.  1  and  2.  As  can  be  seen,  the  agreement  is  good.  The  data 
was  reconstructed  from  the  Fourier  coefficient  given  by  St.  Hilaire  and  Carta 
(NASA  CR-165927).  The  third  measured  data  point  on  the  pressure  surface 
(x/c  *  .066)  gave  very  erratic  results  and  was  not  plotted  for  most  of  the 
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comparisons*  The  excellent  comparisons  shown  in  Figs*  1  and  2  give  evidence 
to  the  time-accurate  calculation  for  the  surface  pressure* 

Figure  3  presents  a  comparison  at  a  ■  17.7°,  a  >  0.  This  is  near  the 
incidence  where  stall  would  first  be  inferred  from  the  lift  and  moment 
curves.  The  figure  shows  some  discrepancy  between  predicted  and  measured 
values  as  the  data  presents  some  evidence  of  a  vortex  being  shed  on  the 
suction  surface  leading  edge.  The  discrepancy  increases  in  Fig.  4  where  the 
data  clearly  indicates  stall*  The  calculated  plateau  on  the  suction  surface, 
x/c  m  .15,  seems  to  indicate  a  vortex  being  initiated.  Furthermore,  the 
calculated  maximum  suction  peak  at  a  ■  19. 5 8,  Fig.  4,  is  considerably  less 
than  that  at  a  ■  17.7°,  Fig.  3.  Based  upon  the  plateau  and  the  drop  in 
suction  peak  the  calculated  distribution  at  19.5°  appears  to  be  beginning  the 
stall  process.  The  data  at  19.5°  is  presented  with  the  calculation  at 
o  “  19.9°,  a  >  0  in  Fig.  5.  Although  these  are  at  different  values  of  a, 
they  represent  pressure  distributions  at  approximately  the  same  incremental 
time  after  stall  is  initiated;  the  distributions  are  remarkably  similar. 
Comparisons  over  the  downstroke  are  given  in  Figs.  6-8.  Obviously,  the  basic 
trends  are  in  agreement  as  a  strong  qualitative  comparison  is  shown  between 
the  calculation  and  the  measured  data.  The  results  of  this  calculation 
particularly  prior  to  stall,  Figs.  1-3,  indicate  the  time-accuracy  of  the 
procedure  at  least  as  far  as  surface  pressure  is  concerned. 

The  second  case  considered  was  performed  under  the  present  effort.  The 
case  was  a  25.5 Z  thick  Joukowskl  airfoil  entering  a  gust  in  which  the 
vertical  velocity  is  0.25  times  the  freestream  velocity.  This  case  was  also 
calculated  by  Gleslng  via  an  invlscid  analysis  (Ref.  14)  and  although  viscous 
effects  are  present  in  the  present  analysis,  it  is  expected  that  the 
calculated  pressure  gradients  should  be  in  reasonable  agreement  if  little  or 
no  separation  occurs.  Calculated  surface  pressure  distributions  from  both 
techniques  are  presented  in  Figs.  9-12  where  T  indicates  the  location  of  the 
gust  relative  to  the  airfoil  leading  edge.  As  can  be  seen,  the  two 
calculations  are  in  good  agreement.  The  major  discrepancies  are  the 
appearance  of  surface  pressure  discontinuities  in  the  inviscid  solution  which 
are  not  found  in  the  viscous  solution  such  as  at  x/c  ■  .25  in  Fig.  10. 
However,  this  is  the  expected  result  of  viscous  effects.  Contour  plots  of 
pressure  coefficient  are  presented  in  Figs.  13  and  14,  where  it  can  be  seen 
that  initially,  the  pressure  coefficient  is  nearly  symmetric,  but  by  T  -  -.10 
it  has  become  distinctly  asymmetric.  It  should  be  noted  that  at  T  ■  -.10  the 


center  of  the  gust  is  at  x/c  “  -.10,  however,  due  to  viscous  effects  the  gust 
front  is  not  a  sharp  line  but  diffuses  over  a  finite  region.  Therefore,  some 
effect  is  felt  at  T  -  -.10.  Velocity  vector  plots  are  given  in  Figs.  15 
and  16.  This  comparison  along  with  the  dynamic  stall  calculations  of  Ref.  11 
serve  to  confirm  the  time-accuracy  capability  of  the  cascade  code. 

Boundary  Conditions 

Previous  experience  in  solving  Navler-Stokes  equations  has  Indicated  the 
Important  role  boundary  donditlons  play  in  obtaining  accurate  solutions  and 
rapid  numerical  convergence.  Improper  specification  and/or  implementation  of 
boundary  conditions  can  lead  to  problems  in  achieving  convergence  and 
solution  accuracy  for  the  solution  procedure.  In  order  to  set  a  framework 
for  the  discussion  to  follow,  it  is  useful  to  re-lew  the  manner  in  which 
boundary  conditions  are  sat  for  calculations  in  which  a  steady  state  solution 
is  sought.  In  these  steady  state  solution  cases,  the  Navier-Stokes  cascade 
analysis  followed  the  suggestion  of  Briley  and  McDonald  (Ref.  15)  which 
specifies  upstream  total  pressure  and  downstream  static  pressure  conditions. 
For  the  cascade  system  shown  in  Fig.  17,  AB  and  CD  are  periodic  boundaries 
and  periodic  conditions  are  set  here. 

Specification  of  upstream  and  downstream  conditions  is  somewhat  more 
difficult.  For  an  isolated  cascade,  boundary  conditions  for  the  differential 
equations  may  be  known  at  both  upstream  infinity  and  downstream  infinity. 
However,  since  computation  economics  argues  for  placing  grid  points  in  the 
vicinity  of  the  cascade  and  minimizing  the  number  of  grid  points  far  from  the 
cascade,  the  upstream  and  downstream  computational  boundaries  should  be  set 
as  close  to  the  cascade  as  is  practical.  However,  when  the  upsteam  boundary 
is  placed  close  to  the  cascade  upstream  influence  of  the  cascade  modifies  the 
flow  such  that  far  field  function  boundary  conditions  may  be  inappropriate. 

In  the  approach  used  to  date,  the  suggestion  of  Ref.  15  is  followed 
which  sets  total  pressure  on  boundary  BC  (see  Fig.  17).  Unless  boundary  BC 
is  very  far  upstream,  the  flew  velocity  along  this  boundary  will  not  be  equal 
to  the  velocty  at  upstream  infinity  since  some  inviscid  deceleration  will 
have  occurred.  However,  as  long  as  the  boundary  is  upstream  of  the  region  of 
any  significant  viscous  or  shock  phenomena,  the  total  pressure  on  this 
boundary  will  be  equal  to  the  total  pressure  at  upstream  infinity.  Hence, 


total  pressure  is  an  appropriate  boundary  condition  realistically  modeling 
the  desired  flow  condition.  In  addition  to  specifying  upstream  total 
pressure,  it  is  necessary  to  specify  the  inlet  flow  angle.  In  the  present 
calculation,  a  specified  constant  angle  was  assumed  on  the  upstream 
boundary.  The  upstream  boundary  condition  specification  was  completed  by 
setting  the  first  derivative  of  density  to  zero.  The  downstream  boundary  was 
treated  by  setting  a  constant  static  pressure  as  a  boundary  condition,  and  by 
setting  second  derivatives  of  both  velocity  components  equal  to  zero  at  this 
location.  In  the  present  application,  a  constant  static  pressure  was  set  at 
downstream  infinity.  Implicit  in  this  assumption  is  that  the  downstream 
boundary  is  located  in  a  region  where  pressure  is  uniform  in  the  transverse 
direction  although  a  nonuniform  specification  is  permitted  numerically  and 
possible  physically.  The  final  boundary  conditions  to  be  considered  are  the 
conditions  along  the  blade  surface.  Here  no-slip  and  no  through-flow 
conditions  were  applied  leading  to  a  specification  of  zero  velocity  on  the 
surface.  An  inviscid  transverse  momentum  equation  was  applied  on  the  surface 
leading  to  a  boundary  condition  approximation  of  zero  transverse  pressure 
gradient  being  applied. 

The  second  item  which  must  be  considered  in  regard  to  boundary 
conditions  is  their  implementation.  Both  the  upstream  and  downstream 
boundaries  have  boundary  conditions  associated  with  them  which  are  nonlinear 
functions  of  the  dependent  variables.  These  are  the  specifications  of  total 
pressure  on  the  upstream  boundary  and  static  pressure  on  the  downstream 
boundary.  These  nonlinear  boundary  conditions  are  linearized  in  the  same 
manner  as  the  governing  equations,  via  a  Taylor  expansion  of  the  dependent 
variables  in  time,  and  then  solved  Implicitly  along  with  the  interior  point 
equations.  Although  points  on  the  periodic  lines  and  the  branch  cut  are 
boundaries  of  the  computational  regime,  they  are  interior  flow  fields  points 
and  must  be  treated  as  such.  The  present  technique  replaces  derivatives  at 
these  points  by  central  differences.  In  addition,  in  regard  to  the  periodic 
lines  the  procedure  Inverts  a  matrix  with  strict  periodic  boundary 
conditions;  i.e.,  the  periodic  line  values  are  obtained  from  the  implicit 
solution,  rather  than  from  an  extrapolation  or  averaging  procedure  which  uses 
Interior  computational  grid  point  values. 

This  specification  has  worked  very  well  when  steady  solutions  are  sought 
as  convergence  could  be  obtained  very  quickly.  Using  convergence 
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acceleration  techniques  described  in  Refs.  17  and  18  has  allowed  cascade 
calculations  to  converge  very  rapidly  even  for  low  Mach  number  subsonic  flows 
which  in  general  can  be  difficult  to  converge.  A  typical  convergence  history 
Is  presented  in  Fig.  18  which  presents  residual  versus  time  step  number  for 
the  subsonic  cascade  calculation  corresponding  to  the  experiment  of  Hobbs, 
et  al  (Ref.  16).  As  can  be  seen,  the  maximum  residual  in  the  flow  field 
drops  by  three  and  one-half  orders  of  magnitude  in  seventy  time  steps.  The 
reason  for  the  sudden  jump  at  time  sixty  will  he  discussed  shortly.  The 
maximum  residual  is  defined  by  the  maximum  imbalance  of  any  equation  at  any 
point  when  the  time-derivative  terms  are  omitted.  At  seventy  time  steps  the 
solution  has  essentially  stopped  changing  and  as  shown  in  Fig.  19,  the 
surface  pressure  calculated  is  in  very  good  agreement  with  that  measured. 

The  reason  for  the  jump  in  residual  at  time  step  60  is  the  run  protocol.  The 
calculation  was  initiated  with  a  moderate  artificial  dissipation  which  was 
dropped  by  an  order  of  magnitude  at  time  step  60. 

Although  these  boundary  conditions  have  proven  effective  and  accurate  in 
obtaining  steady  solutions,  further  assessment  should  be  made  when 
time-accurate  solutions  are  sought.  Application  of  boundary  conditions  has 
been  an  item  of  recent  high  Interest  (e.g.  Ref.  19).  These  efforts  have 
focused  upon  considerations  such  as  the  effect  of  boundary  conditions  upon 
convergence,  proper  boundary  conditions  for  steady  state  solutions  and  proper 
boundary  conditions  for  time-accurate  solution.  When  the  latter  is  sought, 
boundary  conditions  should  not  impose  nonphysical  wave  reflections. 

Approaches  focusing  upon  nonref lective  boundary  conditions  have  been 
formulated  by  several  authors.  For  example,  Bayliss  and  Turkel  (Ref.  20) 
investigated  boundary  conditions  which  simulate  correct  radiation  of  energy 
outside  of  the  computational  domain.  This  approach  requires  the  boundary 
condition  to  be  a  partial  differential  equation  relating  pressure,  density 
and  velocity.  Applications  for  three-dimensional  spherical  radiation, 
two-dimensional  radiation  and  one-dimensional  radiation  were  considered. 

Rudy  and  Strikwerda  (Ref.  21)  also  considered  nonref lective  boundary 
conditions  and  developed  a  specification  which  for  their  problem  and 
numerical  method  hastened  convergence  to  steady  state,  but  did  not  attempt  to 
maintain  transient  accuracy.  Hedstrom  (Ref.  22)  considered  nonref lective 
boundary  conditions  which  maintained  transient  accuracy.  Other  approaches 
are  given  in  Ref.  19.  Based  upon  the  discussions  in  Refs.  20-22,  a 
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one-dimensional  nonref lective  boundary  condition  was  added  to  the  cascade 
code.  This  condition  relates  pressure  and  velocity  by  the  equation 

3p  -  pc  ■  0 

at  at 

where  c  is  the  speed  of  sound.  The  boundary  condition  was  used  in  a 
calculation  of  the  cascade  of  Ref.  16  and  convergence  was  obtained  although 
no  significant  improvement  was  found  over  the  previous  condition.  However, 
this  may  be  a  more  advantageous  condition  in  transient  flow  situations. 

The  final  item  considered  under  the  present  year's  effort  focused  upon 
choice  of  a  stator-rotor  grid.  Two  possibilities  are  immediately  obvious. 

The  first  would  use  two  distinct  grids  which  overlap.  One  grid  would  be  for 
the  stator  and  one  for  the  rotor.  Each  grid  would  be  a  'C'  or  an  'O'  grid 
such  as  those  currently  used  at  SRA.  Under  this  procedure,  the  solution  for 
the  flow  field  about  the  stator  would  be  marched  through  a  time-step  with  the 
(matching  plane)  downstream  boundary  condition  being  set  by  the  static 
pressure  calculated  for  the  rotor  blade  flow  field  at  the  required  physical 
locations.  Thereafter  the  solution  for  the  rotor  flow  field  would  be  marched 
one  time  step  with  the  upstream  boundary  condition  being  set  by  the  stator 
flow  field  total  pressure  results  at  the  required  location.  Obviously,  such 
a  procedure  requires  the  solution  from  one  region  to  provide  appropriate 
boundary  conditions  for  the  other.  The  major  disadvantage  associated  with 
such  an  approach  is  the  required  lagging  of  boundary  conditions  and  the 
interpolation  to  obtain  boundary  conditions.  The  interpolation  as  well  as 
the  lagged  boundary  conditions  may  adversely  affect  numerical  stability  and 
would  reduce  convergence  rate  even  in  stable  calculations.  Nevertheless, 
this  may  be  a  viable  technique  for  the  stator-rotor  problem. 

The  second  possible  approach  would  use  a  single  grid  in  which  the  stator 
would  be  embedded  within  a  'C'  grid  and  the  rotor  within  an  'H'  grid.  A 
sketch  using  a  geometric  configuration  similar  to  that  of  Dring.  et  al 
(Ref.  5)  is  shown  in  Fig.  20.  With  this  approach  both  rotor  and  stator  flow 
fields  are  calculated  as  a  unit  with  no  need  to  Interpolate  or  lag  boundary 
conditions.  Two  problems  arise  with  this  approach.  First,  the  rotor  being 
embedded  within  an  'H'  grid  may  exhibit  poor  behavior  at  the  leading  edge 
coordinate  singularity,  point  A  of  Fig.  20.  Secondly,  as  the  gap  between 


blades  Is  decreased  the  coordinate  system  will  become  increasingly 
nonor thogonal . 

The  first  of  these  problems,  the  treatment  of  coordinate  singular  point, 
was  investigated  during  the  present  reporting  period  through  a  sample 
calculation.  In  this  calculation,  flow  about  a  Joukowski  airfoil  was 
considered  with  the  computational  grid  being  formed  by  the  potential  flow 
streamlines  and  potential  lines.  This  grid  is  shown  in  Fig.  21.  As  can  be 
seen,  the  grid  is  of  an  'H'  type,  and  has  high  resolution  both  in  the  leading 
edge  region  and  in  the  region  in  the  Immediate  vicinity  of  the  airfoil. 

The  calculation  was  initiated  from  a  uniform  flow  with  an  applied 
no-slip  condition  at  the  airfoil  boundary.  For  the  purpose  of  this 
calculation,  the  flow  was  assumed  laminar  with  a  Reynolds  number  based  upon  a 
chord  of  1.5  x  105,  and  a  free  stream  Mach  number  of  .087.  Convergence  was 
obtained  with  no  problem  and  no  flow  anomalies  were  noted  in  the  vicinity  of 
the  leading  edge.  A  comparison  between  calculated  surface  pressure 
distribution  and  that  obtained  analytically  from  potential  flow  theory  is 
shown  in  Fig.  22.  As  can  be  seen,  the  agreement  in  the  leading  edge  region 
is  very  good.  The  viscous  results  give  a  suction  peak  slightly  less  than  the 
inviscid  results.  The  major  discrepancy  is  near  midchord  where  separation 
occurs  in  the  viscous  laminar,  calculations.  This  obviously  changes  the 
pressure  distribution  from  that  obtained  via  the  potential  flow  analysis  and 
the  difference  is  as  expected.  Velocity  vector  plots  are  shown  in  Fig.  23, 
and  static  pressure  contours  in  Fig.  24.  The  results  obtained  clarly 
indicate  the  ability  of  the  present  Navier-Stokes  code  to  yield  accurate 
results  in  the  vicinity  of  the  'H’-grid  coordinate  singularity.  Therefore, 
it  appears  that  the  ’H’-grid  can  be  used  for  viscous  calculations  of  the  type 
anticipated. 

The  second  problem  concerns  the  gap  spacing.  According  to  Dring,  et  al 
(Ref.  5)  axial  rows  in  turbines  contain  gap  spacing  of  1/4  to  1/2  axial 
chord.  Although  conclusions  must  be  tentative  until  demonstrated,  it  appears 
that  this  gap  spacing  range  would  not  preclude  the  single  grid  approach. 
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